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ABSTRACT 


The objective of this study was to examine the relation- 
ship between forward ship motion and the effects of aero- 
dynamic and hydrodynamic disturbances on motions in the 
horizontal plane. It includes course control and stability 
analysis for the unsteered and steered cases (manual and 
automatic) using several sets of operating conditions. 

Hydrodynamic and aerodynamic effects are considered to 
be functions of hull motion, rudder (steered cases), pro- 
peller and wind effects ratio. A dimensionless mathematical 
model has been developed and solved with respect to ship 


axes. 
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Је INTRODUCTION 


In the last twenty years many mathematical derivations 
have been performed in the area of motion of vehicles ina 
fluid medium. Some of these mathematical models have been 
derived using the six degree of freedom equations of motion. 

In this study of six degree freedom equations are re- 
duced to three degrees of freedom with several assumptions. 
Low ship speed is assumed and coupling effects of roll, 
pitch ~ heave motion on to the horizontal plane motion 
are assumed negligible. It is also assumed that atmospheric 
motion is uniform. 

Aspects of the problem are simultaneous consideration of 
hydrodynamic and aerodynamic effects as well as nonlinearity 
in the external disturbances (especially in aerodynamic terms). 
Equilibrium conditions and perturbation equations are derived 
from the mathematical model. Nonlinear parts of the mathe- 
matical model are linearized with respect to equilibrium 
conditions. 

The effects of wind disturbances on ship motions are 
studied. A digital computer (IBM model 360/67) was used to 
solve the transient response in yaw and sway for three dif- 
ferent wind directions (head, beam, stern) and constant wind 
velocities with constant ship speed. Also the Stiga value 
theorem was used to study ship stability. All these com- 
putations were made for the unsteered (zero degree rudder 
deflection) and steered case (with manual as well as auto- 


matic control). 
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When wind velocities are several times the ship speed, 
even moderate winds create difficulty in controlling the 
ship course and stability (at low ship speeds). These ef- 
fects may or may not be controllable, since the hull hydro- 
dynamic reactions are reduced and ship motion is greatly 
effected by wind force, propeller thrust and coupled rudder 
force. This situation can cause unusually large hull drift 
angles extending to between plus or minus 180 degrees. The 
usual hull hydrodynamic force representation is not applica- 


ble beyond 15 to 20 degrees. 
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ІІ. "THE MATHEMATICAL MODEL 


A. COORDINATE SYSTEMS AND SHIP MOTION 

The coordinate system is defined with respect to the 
ship itself. It's coordinate origin ls chosen to be the 
center of gravity of the ship as shown in Figure 1. 

Positive X axis (lengitudinal axis) is directed to the 
forward direction along the center line of the ship. 

Positive Y axis (transverse axis) is directed to star- 
board, perpendicular to the plane of symmetry. 

Positive Z axis is directed downward, and lies in the 
vertical symmetry plane, perpendicular to the X-Y plane. 

Yaw ls a retatienal motion of the ship about the Z axis. 
(Assume that the ship is moving in the positive X direction 
and that a moment is applied about the Z axis tending to 
rotate or turn the ship in clockwise direction as indicated 
in Figure 2 which shows the positive уам). 

1. BEguatsons of Motion (on the horizontal plane) 

The dynamic response in six degrees freedom for the 


motion of a rigid body [7] is given by equation (1). 


x = MLÜ-RV+OW-X, (8402) +Y (PQ-R) 42, (PR4Q) ] 
Ү = m [U- PWeRU*X . (&«PW) -Y . (P^eR^) ez (RQ-B) ] 
2 = M[W-QU+PV+X. (PR-Ó) +Y (B+OR) -Z, (Q°+4P*) ] 


L = оо mee aoe Veen) | 


= 
|| 


QIy+(1,71,)PR+M[Z. (U-RV=QW) -Xc (W-QU*PV) ] 


2 
!! 


RI„+(I,-I,)PQ+M[X, (V-PW+RU)-Y, (U-RV+QW) ] (1) 
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These equations represent the reaction of the rigid 
body. They do not include external forces such as rudder 
force, wind force, aerodynamic and hydrodynamic force of 
water etc. Applied external hydrodynamic and aerodynamic 
forces are considered to be functions of hull motion, rudder 
angle, propeller and wind. 

Assuming that roll, pitch, and heave motions may be 


neglected on the horizontal plane and that x 2. аге 


pr 
zero for. the center of gravity at the origin of the system 


and noting that on the horizontal plane U = iu 4 iv (linear 


velocity), then equation (l) becomes 


Х = mlu - RVI (Surge) 

Y = m[V + RU] (Sway) 

2 = 0 

L = 0 

M = 0 

"E RI, (Yaw) (2) 


These equations are nonlinear. Analytic solution and 
computer simulation are very difficult. But all nonlinear 
terms might be linearized by using the Taylor series expansion. 


Briefly, small perturbations aregconsidqerea for: 


Х = m (u - гу) 
Ү = m (v + ru) 
№ = кі (3) 
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The external force and moments to be included are: [4] 


X = X(hull) + X(rudder and prop.) + X(wind) 
Y = Y(hull) + Y(rudder and prop.) + Y(wind) 
N = 


N(hull) + N(rudder and prop.) + N(wind) (4) 


So N, X, Y represent total hydrodynamic and aerodynamic forces 
and moments which are generated by the ship motion, rudder, 
propeller and wind. The following assumptions permit elimina- 
tion of relatively unimportant terms in the hydrodynamic and 
aerodynamic expressions: [1] 

1. Ship asymmetry with respect to the centerline plane 
is minor (this permits elimination of odd ordered terms for 
the X force expansion and even ordered terms from the Y and 
N expansion). 

2. All terms higher than third order are unimportant. 

3. Rudder force and moment derivatives of higher than 
first order, and effects of rudder angular rate are negligible. 

4. The lengitudinal component of ship speed is much 
greater than the lateral eonponene (u>>v). 

5. Based on experimental and analytical results, non- 
linear hydrodynamic terms in the yaw and sway equation for the 
stability analysis are considered minor, except for certain 
nonlinear hydrodynamic terms in the surge equations such as 
- 

vr 
6. Atmospheric motion is uniform. 


Thus equation (4) may be written as: [1] 


= 5 2 2.2 
Х - суч+С„и +C,vrtC,u 6 pM Б 
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Ү = ΕΙ FH N Us 
Е s 2 
N = A)r+A,uvt+A,urt+A,v SEN uM NOM UL. (5) 


All hydrodynamic and aerodynamic constants are defined in 
the appendis. l 

The hydrodynamic effects, forces are defined as a 
function of the acceleration and velocities of the hull. 

The rudder deflection will also produce an additional 
resistance, an axial force proportional to =" This force is 
expected to be small. The lateral force and moment are pro- 
psetional to S. [5] 

Forces and moments due to wind acting on the above 
water protions of the hull act on the ship in the horizontal 
plane. It is necessary to determine the relative wind speed 
and direction with respect to the ship in it's own reference 
frame. 

Magnitude and direction of wind is shown on Figure 1. 


The relative wind speed components along the ship axes are 


ua = оу, cos (pity) in X direction 
Ya vel. sin (pity) in XY direction 
Then 
L. 2 E - -1Уа 
U, = (u, + Уз) = Уа = tan — (6) 
Ma 


U. and y, are wind velocity and direction along the earth 


axis. 
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Origin of the > Y 
Ship referente frame 


Oriain p the 
Am | n Fixed. Frame. 


Figure 1. Relation of Wind Velocity and Direction to the 
Ship Axes. 
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From equations (3), (4), and (5), accelerations in 
yaw, sway, and surge are determined: 
So 

X = u(m-xu) - юру 

Y = v(m-yv) + mru 

N = г (1, = Nr) (7) 

The relationship between velocities of the ship with 
respect to the earth axis and velocities of the ship with 


respect to the body axis is: [3] 


Y (t) = u sin (y) + v cos (y) 
Х (t) = u COS (y) - v sin (y) 
also 
y(t) =x (8) 


Where y is the heading angle of the ship relative to the 


original path. 


ог! ginal path 


e. o A = 
> > о А о = -— —— A—— Өү л” wa w Аи 





1] без path 


Figure 2. Positive Direction of Yaw. 
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О 


Gee ee А. сөс ())) ас - Y (0) 


|| 


Х 


i fucos Ου sin (v)) d& + X (о) (9) 


when y is small: 


сип), cos (y)ell О=Џе (equilibrium case) (10) 
So 
У (Е) = бе р + у 
X (t) = Ue -ν.ψ C i 


Combining equation (5) and (7) and putting in nondimensional 


form. [2] 


р.з patata Dira ' ı 22 2. ' ' ta 1112 1112 
5^ (m'-xu')W 5% [ (x FX а а ó^)u +, (Хе +m аа ] 
Po 3 ШС б бар ! iS i ! 2 12-112 
> (m v ) V 2% [y ¡UV+2 (y, m')ur+y'öu ΝΟ] 


Pius ' . eee. 3 ' ! ' 2 'N T2 
e a 9 == 4 
5 (I, No ΤῈ 2* [N uv AN ur+N du te NUS (12) 


where 


Га, 3g2N'zN ae, 32%! 
2 аа а 2 а 


NE М! 
Tpm ! ' 2 a) У а2 5 
N^ DUM (> sin(2y_)+, sin(y.)) 
Pa 252! Do 2712 
NEL ap m у 
(71172 V 5 
y- 029: ann) 


да 


g2yu2x'=xK =Pe2y2x! 
2 aa a2 a 


A ZN 
K3 SAM сео) (ЧЕЗ) 
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' = = l >> 

P^ ZU ; usU since u>>v) 

Expanding in Taylor's series around the equilibrium conditions: 
u = ue v = ve é = бе 
F=x=x=u=ve= =0 (14) 


Substitute equation (14) into equation (12) 


N'ueve + N' ée ue? + N'ue? = 0 
M Ó a 

! + ! 2 + ! 2 = 
y„ueve y ,9e ue y, ue 0 


tue? + ! 5252 2 + s! 2 + x'ue? = 
x ue Хх; óe“ ue ха це каше 0 (15) 


Sideslip velocity and rudder angle in the equilibrium condi- 


tion can be expressed as: [1] 


Y 'N.' =  'Y.' 
a 


ve! = —2 6 Š 
М, τς - Y, N. 
МҮ τ Υ "№ ! 

бае = ToU — — И (16) 
'  _ ' 1 
Ny ш; ἵν но 


Perturbation terms (small deviation from the initial straight 


course) may be denoted by bars. So actual motion is 


ue = ue + U V = ve + V с = бе + S ф = 


е | 


0 = 0 ú = ú у= ү. (17) 


Perturbations in aerodynamic quantities due to Y, v' and u' 


ares 
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Where 


oN! 


9} 


oN! 


av! 


oN! 


ді! 


әү’ 
9х 


o ' 


ον ' 


ON 


au! 


_ wo au! 

N! = — ç + — 
a t 
дф ду 

4 Qu. ~ οὓς 

NE Ξ __а ф + кые. 
а дү dy! 
ox! ох. 

дү) 9 V ' 


= ΡΞ 294 


-2Uu' sj ΜΕ. ' 
20. sin Y, 5 sin 2%. + N 


+ N' 
а2 


+ 


oN ) 
ді! 
дү! 

а 
911! 
ox! 

a 


(18) 
au! 


t2 1 
ЕР cos ψ ο) (ЈА чл) 


а, sin Y.) 


М! 


' ' ' _ај : 
ΠΝ, cos 2. + поз ψ ο) + -— sin 2. 


sin V ο) sin Ya! 


= ! ! ! ! -cj 
рд Utm cos 2% а + ЧА cos Фа) ( біп ψ ο) 


+ 2(—1 sin 20, + Ni, sin Vi) cos y] 


2 


Es ' ' v2 a = t аз 2 
ρ Ya, [cos De UM us) 205 біп Фа] 


а 


m crgo σι en) 
a, UL E (1 + sin Va? 


= р! 


' ' 
a БЕ =, 


С біп 
os ұға un 





px 


а _ _ ! ! : ' 2 ' 
ax' 
a ' ' ' : 
5 -ρ! Σ' U! cos uw y 
av! а а] а а а 
IX Е 
MES = ! ! ! 
E. pP} x, Us (1 + cos Ú.) (19) 


The coefficients in equation (19) where taken from reference 
[2]. Aerodynamic coefficients are shown in Figure 3. 
Substitute equation (17) and (18) into equation (12) 


and obtain the perturbation equation in the following form: 


. _ 2х. _ 2х. _ к 2X ^ 
жес. ч b= ! ! lice и І ! ' 
(m Xe )u (oe +m IN P d 47 +(2X +2X5 ¿0 + | ) 
V 2u 
u Τ2Χςρ ὃ Š 
Ñ 2 2 ay! 
n. ' i Em 4 i 1 ' CS ' 
I A A Y Li. О θες) 
x 2v 2u 
џудо 
τ - 283. 2м. _ 2N! 
(I!-N:')r'zN'r'*——Jc-(N'*——)v'c-(N'v'*2N!ó6 4——) 
Z y ¥ 2 V N. ve ó e y! 
ч'+М 6 ` (20) 


The mathematical model which is shown in the above equations 
is developed only for the forward ship motion. The stern 
ship motion equations are written by putting negative signs 


in front of some of the coefficients. 
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Exgure 3, Hydrodynamic Coefficients Vs Rudder Force Coefficient. 
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Example: N, (derivative of hydrodynamic yaw moment with 
respect to sideslip velocity) takes a negative value for the 
forward motion, positive value for the stern motion. 

Y, (derivative of hydrodynamic force component in 
Y-axis, with respect to yaw rate) takes negative value for 
the forward motion, positive for the stern motion but is not 
a sensitive parameter since m' is the predominant term in 


E - m). 
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III. СОТКО STUDIES 


The linearized equation of motion was derived and found 
in equation (20). In this part two different problems were 
simulated. One of them was manually steered ship motion in 
wind disturbance. The other was automatically steered ship 
motion under the same effect. So computer program (II) and 
(III) were developed for these purposes respectively. Тһе 
computer model is shown in Figure 4, for these programs. 
Program (II) and (III) are similar to each other. There is 
only one difference in section 4 of both programs. Program 
(II) was integrated using the second order equation, however 
program (III) was integrated with the third order equation. 
This difference was caused by the feedback loop (page 58 ). 
Also the Laplace Transform Technique could have been used to 
obtain the system description for simplicity and system 
characteristic equation for the parameter plane program. The 
characteristic equation was calculated in the computer pro- 
gram. Also three more additional and different mathematical 
sections were calcualted in the same program. One of these 
sections included the parameter calculation from equation 
(6) through equation (20). Secondly the transfer functions 
Of yaw and sway were calculated from the three degree of 
freedom equation. The last section was the polynomial cal- 
culation for the eigen values. 

The last program, (parameter plane (IV), was used to get 
optimal values for yaw gain constant and yaw rate gain con- 


stant using the constant zeta (function omega) and constant . 
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— (function zeta) curves for automatically steered ship 
motion. According to reference [1] and [6], these optimal 
values were chosen from Figures 40 through 44 assuming 
zeta = 6, omega = 11.5 for varying wind direction and 
velocities. 

The transform of the linearized equations of motion (20) 
from the T domain to the S domain, for all initial conditions 


equal to zero are: 


(ABC. 82-AAD. $-BFA)Y- (BGD.8) -(BGC. 8) IAS 
6 


© |< 


(-86Ε.6-ΒΕΡ) ᾖ- (Ώ6ς.5) 5 +(АВЬ.52-ВСЕ.5) = ААС. 5 (21) 
ó 5 


(-CAA.S-BFG)U- (ABE. S2-CAB.8)9. - (BGA.5) - » CAC.6 
5 5 


EemuBtron (21) different letters are used for simplicity. 


' 
For example: АВС = І! - Ne' , AAD = Ni, BFA = 2a , 
АВЕ = m' - хы! p etec. 
let 
an C ; E B 
S S 
Applying Cramer's rule, the transfer function for И and В 
5 5 
аге: 
ААЕ “BGD. S AAE 
4 ААС -BGG.S AAC 
Y 
Е | САС ABE.S?-CAB.S САС (22) 
$ (АВС.52-ААр.5-ВҒА) -BGD.S -BGC.S 
(-BGE.S-BFD) -BGG.S ABD.S?-BGF.S 
(-CAA.S-BFG) ABE.S?-CAB.S  -BGA.s ‚| 
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The open loop transfer function is: 





» VVA.S*+VVC.S+VVF 

ó SSA.S"+SSB.S3+SSE.S°+SBS.S+SES (22) 
Where 

VVA = - ABC.ABD.ABE etc. 


All terms were calculated using program l. 
In this case, the results are a second order numerator 


and fourth order denominator. 





Transfer function БАӘ), was also obtained using Cramer's 
5 ($) 
rule. 
ABC.S2-AAD.S-BFA -BGD.S AAE 
-BGE.S-BFD -BGG.S AAC 
B -CAA.S-BFG ABE.S?-CAB.S САС 
== ——— —— RÜ -—— MO (23) 
EI ape ο. Ιμπν —B@p.s -BEC е 
-BGE.S-BFD -BGG.S ABD.S?-BGF.S 
-CAA.S-BFG ABE.S?-CAB.S -BGA.S 


The transfer function between B and $ is: 


B WWA.S?-WWC.S?-WWF.S-WVB 
о (24) 
$ SSA.S°+SSB.S'+SSE.S3+SBS.S7+SES.S 


WWA = -ABC.AAC.ABE etc. 


Again all terms were calculated using program l. 
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For the lateral displacement, sway is defined by 





WWA. S3+WWC , S2+WWF . S+WVB 
SSA.S°+SSB.S'+SSE.S3+SBS.S7+SES.S 


o 


VVA.S?4VVC.S-VVF 
SSA.S"+SSB.S3+SSE.S?+SBS.S+SES 


f 





ie transfer function is 


ү WWA.S +(WWC+VVA)S +(WWF+VVC)S+WVB+VVF 


— 
Dam 


S NECS C ULTENP SUTUOSETOSBA)S TUSBSTSSB)S t(SESTSSE)S ἘΒΡΕΘ.5ΠΘΕΘ.8 
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A.  MANUALLY STEERED SHIP MOTION 
As mentioned before, the goal of this section is to get 
the ship trajectories in constant wind effect with low ship 
speed. For this purpose computer program (II) was developed 
and applied to equation (20). Two situations were considered: 
Situation 1: zero degree rudder deflection (unsteered). 
Situation 2: positive 15 degree rudder deflection. 
All these situation were studied assuming that the ship moved 
in the forward direction a а 1 mile/hr speed in.a calm sea. 
Each of the wind velocities was applied suddenly. 
Computer program II responses are shown in Figures 5 
through 30 for yaw and sway. These results were also tabu- 


lated in Tables 1 through 4. The first test on the unsteered 
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situation is shown in the computer outputs trem Figure 5 to 
Figure 24 for three different wind directions and various 
velocities. Tabulated values are shown on Table 1 and 2, 
for yaw and sway. 

The effects are examined by referring to the tables. 

Casel: Head Wind 

As is shown in Tables 1 and 2 and Figures 5 - 10, effects 
of the head wind were small on the straight forward ship 
motion. In this test, a 5 mile/hr wind velocity is applied 
to the ship which created ο; rad. variation on the course 
and drift was -7.4х1073 fc. 1015 minutea, 

Case 2: Beam Wind 

Secondly, beam wind disturbances were examined for various 
velocities. According to Tables 1 and 2 and Figures 11 -16, 
the ship greatly deflected from a straight course where the 
wind velocity exceeded t miles/hr. Deviation was positive 
22 degree (.395 rad) on the straight course and positive 
9.2 ft drift in 5 minutes when the 7 miles wind velocity was 
applied, as shown in Figures 15 and 16. 

Case 3: Stern Wind 

In this part, the ship was disturbed with stern wind 
which created great difficulty on the straight course. Five 
miles/hr wind velocity created -2.29 degree (-.04 rad) devia- 
tion and -.28 feet drift in 5 minute. But 7 miles/hr wind 
velocity causes the ship to deviate by -332.3 degree from 
the straight course in 4 minute 10 second. Drift was -110 ft 


in 4 minute but this value was reduced to -15 ft in 5 minute. 
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κοιν the system was unstable and the ship was turned 
around, There is no chance to control the ship. Drift was 
negative 110 ft in 4 minute but this value was reduced the 
negative 15 ft in 5 minute. The summarization of this test 
can be stated as the following: computer program (II) re- 
sponses were tabulated for 3 different wind directions and 
various velocities. According to the table, wind distur- 
bance was not so important on the course from any direction 
until 5 miles/hr wind velocity was reached. Critical points 
were 5 and 7 miles/hr wind velocities for stern and beam 
wind respectively. 

Second test: 15 degree rudder displacement was applied 
under the beam and stern wind disturbances. The main pur- 
pose was to show the effect of the rudder displacement on the 
ship course with wind. As is known, positive rudder displace- 
ment causes negative yaw. All responses are shown in Figures 
25-30 and tabulated values are shown in Tables 3 and 4. Ас- 
cording to the Tables 3 and 4, when the ship is effected 5 
miles beam wind velocity and positive 15 degree rudder dis- 
placement. So, ship was deflected -327 degree from the 
original path in 5 minute. Drift was 80 ft in 2.5 minute 
and -100 ft in 5 minute. As mentioned earlier, 5 miles/hr 
beam wind with zero rudder displacement has created 8 x oe 
rad. deviation from the original path, as shown in Table 1l. 
But the same condition with 15° rudder displacement has -327° 
deviation from the original path. In this case, the ship 
was so effected with rudder force and started to turn around 


the path. 
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Gure WD DI, (degree) WD.VEL. (miles) Time(sec) Yaw(rad) 


180 3 300 -4 
180 5 300 -18 
90 5 300 -5.7. 


Table 3. Variation of Yaw for Manually Steered Ship (S=15°). 


соке WD.DI. (degree) WD.VEL.(miles) Time(sec) Sway(ft) 


85 9 
180 3 {2490 Ü 1g 
180 5 nonlinear oscillation 
150 80 
τμ E 300 [100 


Table 4. Variation of Sway for Manually Steered Ship (S=15°). 
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The other situation is a stern wind disturbance with 15° 
rudder displacement. Also known from Table 1, 5 miles/hr 
stern wind velocity with zero rudder displacement has created 
-.04 rad. deviation from the original path. In this case, 

15 degree rudder displacement created -18 rad. deviation in 
5 minute. Ship was turned 360° around the course in 2 minute 
and drifted with divergent oscillation. 

All these tests were made with ship speed being 1 mile/hr 
in the calm water and values of wind velocities were suddenly 


applied. 


В. AUTOMATICALLY STEERED SHIP 

Earlier tests show that the ship was effected by wind 
and rudder. 

In this section, another test was developed to minimize 
the deviation from the original path. For this reason 
rudder deflection is used and is written as a function of v. 


The equation for the rudder deflection is: 


5 = а 5d + Sd = κ, ф + ky 


where 
a = rudder response time constant (time lag) 
k = yaw gain constant 
k. = yaw rate gain constant 


Thus, the system would be: 
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These control effectors are the function of some parameter 


that could be measured by sensors, usually a motion parameter. 


VVA. S + ММС. 5 + ММ Е 
SSB.s*4,SSb.S*4 SSE.S'4 38 S. SE SES 








Let 
m uo 
a 1 
с = Ë _ __Wa.s?+ VVc.S+VVE m (26) 
6d SSA.S*+SSB.S*+SSE.S?*+SBS.S+SES Sta, 
H = K + K,S 
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So the characteristic equation is 


| SSA. S** (SSB*a., /SSA) . S** (SSEta, SSB-a, .VVA. K, ) S? 


2 
e E 9 9 = ә Ф T ° 
+(SBSta, SSE a, VVC K. A VVA.K)S“+(SES а 585 


-a..VVF.K, -a . VVC. K) $* (a, . SES-a 


1 t 91 .VVF.K) = 0 (27) 


1 


For simplicity, let  B,-SSA, B,-SSBto *SSA etc. 


1 1 


5 4 - 
8,'5"+8,5 +(B,-B 


3 38 Е = 2 
2 KS + (В, B +K В-К)5 


4 6 t 


* (Bg-B . K -Ві0:К)5%(811-В •К) = 0 (28) 


ONE 12 


This characteristic equation was calculated using computer 
program l. Computer program III was used to get transient 
responses in yaw and sway for the automatically steered ship. 
Each curve of yaw or sway was developed with different param- 
eter values (for K and Ki), for different wind speed and 
direction. 

The problem is to minimize the variation on the ship 
course. For this reason, the parameter plane program is ap- | 
plied to the system to determine the best value of yaw gain 
constant and yaw rate gain constant. 

The responses of these are shown from Figure 40 to 
Figure 44 for several situations. The values chosen are 
the minimum values shown on the curves according to refer- 
ence [6]. 

The parameter plane outputs are developed by applying 


equation (27) to program IV. Table 5 shows the best values 
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of the yaw gain constant and yaw rate gain constant for the 
automatically steered ship. 

All responses of the automatically steered ship are 
shown in Figure 31 through 33 and tabulated on Tables 6 and 
m 

According to the Tables (6 and 7), the variations on 
the course is minimized for the low wind velocities. In 
the head wind disturbance, course variations were reduced 
БӘ 0 rad. in 5 minute. 

One mile/hr wind velocity response reached 0 rad. in 
20 sec., 3 mile/hr gave 0 rad. response in 65 sec. But 5 
miles/hr wind velocity response was 1x10 4 rad. in 5 minute. 

For the beam wind effection, 3 miles/hr wind velocity 
was applied to the ship. Deviation was 0 гай. іп 70 sec. 

A stern wind was also tried, 1 mile/hr wind velocity 
caused a deviation which was corrected to 0 rad. in 65 sec. 

As is seen above, the ship trajectory was kept on the 
original path. This desired course or original path was 
assumed to be due North. Also the same test condition was 


applied. 
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Figure WD.DI. (degree) WD.VEL. (miles) 


Table 6. 


Moure 


Table 7. 


180 


90 


Time (sec) 


10 
300 


50 


1300 


{ 15 
300 


20 
(300 


30 
(300 


{ 


{ 


{ 


{ .85x10 


{ 


Yaw(rad) 


-1x10 
“0 


-1x10 
-0 


-1x10 ` 
“0 
=) 


1x1074 
.5Х107 


4 


4 


4 


4 


4 


Variation of Yaw for Automatically Steered Ship. 


WD.DI. (degree) 


180 


0 


90. 


WD.VEL. (miles) 


Time (sec) 


300 
300 
300 


300 


Sway (ft) 


10x10 


13x10. 


1l.5x10 


-5x10 


4 


3 


3 


Variation of Sway for Automatically Steered Ship. 
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ТУ. STABILITY ANALYSIS 


Stability may be defined in the following way. Estab- 
lish an equilibrium situation, disturb the equilibrium con- 
dition with arbitrary disturbances. If the ship returns or 
tends to return to the original straight-line condition of 
equilibrium after a disturbances of the smallest amount 
(infinitesimal disturbance), it is stable. If it departs 
or has the tendency to depart from the original straight- 


line equilibrium condition, it is unstable. 


The unstable ship will go in to a starboard or port 


Unstable 
Ship path 
κ” 
7 


Ship path 
unstable 
Ship path 
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Stability analysis is performed by solving for the eigen 
values of the system matrix (4 x 4) which is defined on page 
46. Reference [3] contains straight forward algebraic solu- 
tion of the linearized equations. According to reference 
[3] a complete algebraic solution of linearized equation (20) 
is obtained with a differential operator. Solution result 


of any independent variable E is defined as: 


git gt G t 
Xj (t) = ES + 5215 + .... ΝΠ 
m git 
X. = 2 C.. e (29) 
Bec) i=] 1) 


Where Ei are arbitrary constants depending upon the initial 
conditions and σ1, 02,..... 0, аге the roots of the deter- 


minant of the coefficients in the linearized equations of 


motion. 
In general Ој, Dor .... mn are different in value, the 
ОДЕ Ot. 
terms Cie Е i с-е | К а-ы» c e?m* can not negate one another. 


The only way that each term can go to zero with increasing 
time is for each of the exponents to be negative. Ifo is 


a complex number in the form o=a+tib, the following relation- 


et datib)t а 
=e =e 


ship hold e С (соз bt+i sin bt) and the condition 


1 τη ° ° °. s 


be negative as they are complex numbers (the imaginary oc, 


for stability requires that the real parts of c 


parts of number indicate the angular frequency of oscillation). 
The Routh oriterian can also be applied to determine system 
stability. If we look at equation (20) again, the automat- 


ically steered ship equation would be: 
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ABC.V-AAD.V-BFA.Y-BGC.V-BGD.U = ARE. 


un 


-ВСЕ.ф-ВЕрф +ABD.V-BGF V-BGG U 


= AAC.S 
-CAA.Y-BFGY- BGA.V+ABE.U-CAB.U = CAC.S 
Q0Q.ytQQO.y = 000.5 + S (30) 
Where 
ООО = к, QOQ = k, QQQ = tr = а 


Equation (30) will be written in the following way. 














-AAD -ВҒА 0 -ВСС 0 -ВСр 0 -AAE 
SESEI. LSBED)  |BMBD|. |-BGE| | O |. I-BGG| О | • |-ААС | _ 
W+ + V+ V+ U+ U+ E* 5-0 
-CAA -BFG О -ВСА АВЕ -САВ 0 -САС 
| 
000 000 0 |4 0 0 0 -000| |-1 
| | (31) 


As is seen from equation (31) there are four independent 
variables namely yaw, sway, surge and the rudder angle. The 


form used in equation (29) is applied. 


m | τ... | μμ. . Р _ м 
Vei eit, ver %,е01%, бы) Time io, Ξ-Σ A efi" (32) 
i=l ^ i71 1=1 1=1 
Where МАГ Vi U; апа A; are arbitrary constants. 


Substitution of equation (32) into equation (31): 


ABCO „RAD “© unn -BGC -BGD -ААЕ 
-BGEO , -BFO _ |һвро,-В6Е| |-всб _ |-ААС _ 
Vut V, + U.* А, 0 
-СААО ; ВЕС -BGA ABEO ‚САВ “(САС 
0000. %0200 0 0 -0000,-1 


(33) 
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Solutions of equation (33) are possible for only values of 


σι which satisfy equation (34). 


(ABCO, “-AADO,-BFA) (-ВбС) (-BGD) (-AAE) 
(~BGEo , ~BFD) (ABDo,-BGF) (-BGG) (-ААС) | 
(-СААс,-ВЕС) (-ВСА) (ABEc;-CAB) (-CAC) 
(0000. +000) 0 0 (-9006.-1) 


(34) 
Program I gives the solution of this equation as fifth 
order polynomial. For different situations such as wind 
velocity and direction other polynomials were also found. 
As is pointed out in page , it is important for the real 
root to be negative for system stability. 
This test was examined in two situations: 


Situation 1: constant wind directions with varied 
velocities. 


Situation 2: constant wind velocities with varied 
direction. 


The first situation was divided into 3 sections: 
Section a: constant head wind with several velocities. 
Section b: constant beam wind with several velocities. 
Section c: constant stern wind with several velocities. 
The second situation was divided into 2 section: 


Section a: constant 7. miles/hr wind velocity with 
several directions. 


Section b: constant 9 miles/hr wind velocity with 
several directions. 


Ша Order to examine the stability of an automatically steered 


ship in wind, all these situations are solved with computer 
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program I, tabulated on Tables 8 - 12 and drawn on the 
CALCOMP plotters, from Figures 45 through 49. 
Situation l: This test was developed, constant wind 
direction disturbed the ship. Wind velocities changed 
each time. 
a) Head ας disturbance: 
As is shown in Table 8 and Figure 45, the ship 
would be stable in varied head wind velocities until 13 
miles/hr wind velocity. After this range, the ship could 
be unstable on the horizontal plane. 
b) Beam wind disturbance: 
In this case, the stability range was shorter 
(from 1 to 8 miles/hr). After the 8 miles/hr wind velocities 
the ship could be strictly unstable, shown on Table 9 and 
Figure 46. 
с) Stern wind disturbance: | 
Stern wind effect on the stability might be the 
shortest range, as shown on Figure 47 and Table 10. This 
range included from 1 to 6 miles/hr. After 6 miles/hr wind 
velocities, the ship could be unstable. 
Situation 2: In this case, constant wind velocity dis- 
turbed the ship. Wind direction was changed for each test. 
a) In 7 miles/hr wind velocity with varied direction. 
From 0 degree through 180 degree wind directions were tried 
in 7 miles/hr wind velocity. According to Table 11 and 
Figure 48, the ship could be unstable when the wind dis- 


turbed the ship from 30 degree, 105 degree and 180 degree 
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wind direction. Іп the other wind direction and 7 miles/hr 
wind velocity the ship.trajectory could be stable. 

b) In 9 miles/hr wind velocity with varied wind 
direction. As is shown on Figure 49 and Table 12, 9 miles/hr 
wind velocity with each direction were applied to the system. 
The ship could be stable when the wind direction was head, 

60 degree - 75 degree and 120 degree - 135 degree. For the 
rest of the directions (0° - 180°) the ship trajectory was 
unstable. 

All of these results show that some of the wind direc- 
tions and velocities played a very active role. These direc- 
tions are 30 degree, 105 degree and 180 degree. After the 
9 miles/hr wind velocities, many of the wind directions had 


adverse affects on stability. 
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Vs DISCUSSION OF RESULTS AND CONCLUSIONS 


The procedures for system course control and stability 
analysis on the horizontal plane were described. Applying 
the datas which are taken from some references, to the 
three degree freedom equation gave us course control and 
stability analysis computer outputs. As a result of these 
outputs some comments can be made. First of all the un- 
steered case was examined. It is seen that the stern wind 
effects the straight course of the ship the most. When the 
wind velocity reaches 7 miles/hr this effect becomes ex- 
tremely high. In 5 minutes, the deviation of the course 
is about 327 degree. At the same time the effect of the 
head wind on the straight course compared to the stern and 
beam wind seems to be much less; which is seen from the 
computer outputs. 

To see the effect of the rudder deflection on the 
straight course there should be positive 15 degree applied 
to the rudder for beam and stern wind. Under these con- 
ditions negative yaw and sway were found. When the ship 
is under stern wind, applying positive 15 degree rudder 
deflection will create extremely negative yaw and sway. 

As it is seen in Figures 29 and 30, when the stern wind 
velocity reaches 5 miles in 50 seconds, 15 degree rudder 
deflection will create minus 145 degree deviation from 


tee course. 
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Computer output results were also obtained for an auto- 
matically steered ship. Parameter plane computer programs 
were applied to the closed loop system for choosing the 
best value of yaw rate gain constant and yaw gain constant. 
According to the automatically steered ship computer out- 
puts, in low wind speeds the course variation is minimized. 
Seay Can't be controlled very well but the value of the 
drift has decreased compared to the unsteered case. 

In‘a final study, using the eigen value theorem, system 
stability and the critical wind velocity are found and 
tabulated.  Computed results show that 30 degree, 105 degree 
and 180 degree wind direction are important on the ship 
course in X miles/hr wind velocities. Also 9 miles/hr wind 
velocity in many direction were so effected. 

All the studies done above are made for a ship which 
has a speed of 1 mile/hr and limited special full loaded 


ship form. 
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APPENDIX 


Hydrodynamic and aerodynamic coefficients where used 


in equation (5) 


= е 5 e l = e 3 • ! = e 3 • ! 
Ај 5 £ Ne! By 5 2 Ye σι 5 2 Xe 

= е 3 ' = е 2 ! = е 4 ! 
А. 5 2 N. B, > ў, Yy с. 2 ў, X 

= е 0% ' = е 3 ! = е 3 ' 
ue ^ ἂν о а TAS 

= е 3 ! ως е 2 ! 26 е 2 ! 
Ag 2 E Ns В, 2 A 75 Сд 2 A Ass 

Where 

А, = Coefficients of hydrodynamic moment in yaw 
в. = Coefficients of hydrodynamic forces in sway 
C; = Coefficients of hydrodynamic forces in surge 
1 = 1,2,3,4 


№, Νο, Yi; X; are coefficients of aerodynamic forces 
and moments 


Aerodynamic forces and moment 
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u 
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A 21212 i ya Ха, 51п фа 
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Figure 5. Variation of Yaw Versus Time for Wind Direction 
0009, Velocity 1 mile/hr. 
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Figure 6. Variation of Sway Versus Time for Wind Direction 
000°, Velocity 1 mile/hr. 
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Variation of Sway Versus Time for Wind Direction 
000°, Velocity 3 miles/hr. 
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Variation of Yaw Versus Time for Wind Direction 
000°, Velocity 5 mlles/hr. 
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Figure 10. Variation of Sway Versus Time for Wind Direction 
000°, Velocity 5 miles/hr. 
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Figure 11. 


Variation of Yaw Versus Time for Wind Direction 
090°, Velocity 3 miles/hr. 
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Variation of Yaw Versus Time for Wind Direction 
090°, Velocity 7 miles/hr. 
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Variation of Yaw Versus Tim 
180°, Velocity 1 mile/hr. 


e for Wind Direction 


70 





Sway (48) 
4.10" 


5 
| E м“ а 
TIME (Sec) 
Figure 18. 


Variation of Sway Versus Time 


for Wind Direction 
180°, Velocity 1 mile/hr. 


71 





0 
с 
с 
nN 
ел 
2 8 
Е 4 
e^ 
> 
3 
| 
νο 
to 
5 
1 
Figure 19. 





т 
о 
fae 


Usu. 
A 


ТІМЕ (Sec) 


Variation of Yaw Versus Time for Wind Direction 
180°, Velocity 3 miles/hr. 
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Figure 20. Variation of Sway Versus Time for Wind Direction 
180°, Velocity 3 miles/hr. 
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Figure 21. Variation of Yaw Versus Time f 
180°, Velocity 5 miles/hr. 
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Figure 22. Variation of Sway Versus Time for Wind Direction 
- 180°, Velocity 5 miles/hr. 
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Variation of Yaw Versus Time for Wind Direction 
090°, Velocity 5 miles/hr, Rudder Deflection 15°. 
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Variation of Sway Versus Time for Wind Direction 
090°, Velocity 5 miles/hr, Rudder Deflection 15°. 
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Variation of Yaw Versus Time for Wind Direction 
180°, Velocity 3 miles/hr, Rudder Deflection 15°. 
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Figure 28. Variation of Sway Versus Time for Wind Direction 
1809, Velocity 3 miles/hr, Rudder Deflection 15%, 
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180°, Velocity 5 miles/hr, Rudder Deflection 15°. 
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Figure 30. Variation of Sway Versus Time for Wind Direction 


180°, Velocity 5 miles/hr, Rudder Deflection 15°. 
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Figure 31. Variation of Yaw Versus Time for Wind Direction 
000°, Velocity 1 mile/hr, Automatically Steered. 
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Variation of Sway Versus Time for Wind Direction 
000°, Velocity 1 mile/hr, Automatically Steered. 
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Variation of Yaw Versus Time for Wind Direction 
000°, Velocity 3 miles/hr, Automatically Steered. 
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Figure 35. Variation of Yaw Versus Time for 
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Velocity 5 miles/hr, Automatically Steered. 
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Variation of Yaw Versus Time for Wind Direction 
090°, Velocity 3 miles/hr, Automatically Steered. 
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Figure 38. Variation of Yaw Versus Time for Wind Direction 
180°, Velocity 1 mile/hr, Automatically Steered. 
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Figure 39, Variation of Sway Versus Time for Wind Direction 
180%, Velocity 1 mile/hr, Automatically Steered. 
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Figure 40. 


93 





Figure 41. 
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Yaw Gain Constant Versus Yaw Rate Gain Constant 
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for Wind Direction 000?, Velocity 5 miles/hr. 
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Yaw Gain Constant Versus Yaw Rate Gain Constant 
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for Wind Direction 180°, Velocity l mile/hr. 
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Figure 48. Real ©. Versus Wind Directton for 7 miles/hr Wind 
velocity. 
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